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Abstract 

A recursive Green's function method is employed to calculate the density-of-states, transmission 
function, and current through a 150 layer (3,3) armchair nanotorus (1800 atoms) with laterally 
attached metallic leads as functions of relative lead angle and magnetic flux. Plateaus in the 
transmissivity through the torus occur over wide ranges of lead placement, accompanied by en- 
hancements in the transmissivity through the torus as magnetic flux normal to the toroidal plane 
is varied. 
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1. INTRODUCTION 



Carbon nanotube structures and graphene materials show promise 

for advancing nanoelectronics in the 21st century. Their unique mechanical, optical and 
electronic properties make them important materials for studying coherent quantum trans- 
port and quantum control of both charge and spin 

Toroidal carbon nanotube structures, also known as carbon nanotori, could play interesting 
roles in nanoelectronics, quantum computing, and as biosensors 19|, l20|]. To date, most of 
the efforts directed towards quantum rings have focused on flat two-dimensional structures 
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M M, \m m m m m lad, lail, laa, m M, m, m, l37|, lasi, l39| . whiie these objects 



can be more easily modelled in comparison to their toroidal counter parts, toroidal struc- 
tures allow for motions around their minor radius in addition to azimuthal motion, which 
could give rise to novel phenomena 40|, |4l|, |42] . Some of the potential applications by which 



the unique transport features of nanoscale tori might be exploited are: 

• as macrosopic, three-dimensional molecular rings with persistent current effects in 
magnetic fields or azimuthal and dipole-type electronic excitations in microwave fields; 

• as molecular Aharonov-Bohm oscillators with modulations in the current as a function 
of angle between the attached metallic leads or magnetic flux parallel to the toroidal 
symmetry axis; 

• as biosensors, in which modifications of device properties are altered after covalently 
attaching a biopolymer. 

In this paper, the density-of-states D{E), transmission function T{E), and current 
through the torus Isd between two attached metallic leads under small bias are calcu- 
lated via tight-binding Green's function methods. The nanotorus modelled here consists of 
150 layers of a (3,3) armchair unit cell giving 1800 carbon atoms. It measures 4A in width 
and has a central diameter of D = 116A, making it thin compared to the graphene ring 
discussed in [43| where the ring width-to-radius is set to 1 : 5 or 1 : 10 respectively. The 
opening angle between the leads is varied along with magnetic flux $ through the torus 
to study oscillations of Isd due to constructive or destructive interference of the electronic 
pathways through the left and right semi-circular toroidal branches. Electron-charging ef- 
fects with band-bending and modification of the Fermi-level in the torus have been shown 
to be small in similar small-bias device architectures utilizing carbon nanotubes, or can be 
sufficiently approximated with a linear potential profile in the device region without major 



change of results compared to the exact self-consistent treatment [4J]. Thus, a full self- 
consistent treatment of electron-electron correlation and Coulomb repulsion effects on the 
quantum transport will not be discussed here. The conductivity G and Isd for the small-bias 
case can be obtained from T{E) by a convolution over E with the difference in the Fermi 



distributions functions fi^R of the left and right metalhc leads [45 



This paper is organized as follows: In section 2 the Green's function formalism in the 
tight-binding approximation by which the results obtained here are presented. In section 
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3 graphs of T{E) and D{E) and of I so as functions of different small-bias source-drain 
voltages VsD are shown. A magnetic field parallel to the central symmetry axis (z-axis) 
is applied to study as a function of magnetic flux through the structure. Section 4 is 
reserved for conclusions and suggestions for further work. 



2. TRANSPORT IN THE TIGHTBINDING APPROXIMATION 

The device setup for the toroidal nanotube is sketched in Figure 1. Each semi- infinite 
metallic leads touches the toroidal surface at four atomic contact points. The positions of 
these contacts can be determined from the toroidal vector r 

y{9, (p) = {R + a cos6)ep + a sinOez (1) 

and the relative positions of the two leads. Effectively, the curvature of the torus is neglected 
at the lead interface. The opening angle a between the two leads can be varied from 180° 
{back-to-back configuration) down to 45° (or 315° respectively). 



The Hamiltonian for a nanotube in the tight-binding approximation is 45 



with Ei the on-site energies and the tij the hopping terms. Here all Ei=0 and the tij for 
the torus will be set equal to v = — 3.1eV in accordance with the literature j46|. These 
assignments neglect any curvature effects on the tij. 

The device is modelled as an open system under a small bias voltage using a Green's 



function method |45l. 1471. 148l. I49l|. The semi-infinite leads can be absorbed into the definition 
of a device Green's function Gd that serves to yield transport properties of the nanotorus. 
Self-energy corrections S describing the coupling of the device region to both leads can be 
folded into a device Hamiltonian Ha that carries the retarded and advanced Greens functions 
gI^''^\e) as solutions to a unit source via 

[E-Hd-J:L-^R± ^V] G?''^ = /. (3) 

The Hamiltonian Hd is an 1800 x 1800 matrix at the C-C level but may be re-written as 
a 150 X 150 matrix whose elements consist of 12 x 12 intra or interlayer armchair coupling 
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submatrices A and V respectively, 
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A numerically fast, recursive algorithm to invert the bracketed term in Eq.(3) to compute 
Gfi is known for an infinite nanotube. For the toroidal carbon nanotorus, however, there is an 
additional complication due to the corner elements in the Hamiltonian needed to achieve ring 
closure. An algorithm that extends the simple case to include toroidal ring closure has been 
developed by one of us [s^ . and is employed here to determine the device Green's function. 
Tjl^ji can be calculated analytically for semi-infinite metallic leads. All relevant transport 
properties can be derived from the retarded and advanced Green's functions G[^''^^ (E) of the 
device region. 

The Greens function of the semi-infinite leads, Qm, can be determined analytically. Dis- 
cretizing along the longitudinal direction with lattice spacing a as is standard practice 25|, |45 
gives for a lead with cross-sectional area LyLz 
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t ma^ 
with Emn the standard particle in a box eigenvalue expression and Ep the Fermi energy of 
the metal, included to account for the zero of energy being set at the tube Fermi energy 
E = Q [45|. Eq.(5) is understood to be evaluated on the layer closest to the torus. 

The exact value of the corresponding electronic coupling of the first atomic layer 
of the metallic lead surface to the torus surface is strongly dependent on the choice of the 
metal and its lattice structure. = —0.25 eV is chosen as initial default [45| and its value 
later modified and its effect on T{E) later studied. 

The formula for D{E), T{E) and Isd under bias are summarized below. For this, the 
coupling terms Ti jilE) for the lead-torus connections are relevant and computed from the 
self energies and subsequently from gm{E): 



T,{E) = 2m - Sl(E)J = 2nV^Im [g^iE)] V^, k = L,R. 

From Eq.(3), (4) and (7), D{E) and transmission T{E) are determined as: 

D{E) = G, (r^ + r^)G^; 



(7) 



(8) 
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T{E) = Trace TLGdTRGl 



(9) 



The integrated source-drain current Isd is finally given as: 



dE T{E) [fo {E - ^l) -fniE- 



(10) 



— oo 



with foiE) 



1 



[exp(£;-/.)+l] • 



Finally, we are also interested in studying coherence in electron transport induced by 
electrons interacting with an external magnetic field. The vector potential A (r) for a static 
magnetic field aligned with the toroidal symmetry axis can be written in the Coulomb gauge 
as: 



It can be shown that A (r) can be incorporated into the in terms of a phase factor ( via 



3. RESULTS 

A. Density-of-states D{E) and transmission function T{E) 

In Figure 2, density-of-states D{E) and transmission function T{E) plots are depicted 
for different choices of Bq and tj^Qp- In Figures 2a and 2b, results for D{E) and T{E) are 
compared at Bq = 0, Bq = 0.5 T and Bq = 1.5 T. In each case the overall peak structures of 
D{E) and T{E) stay essentially unchanged. Van Hove peaks are clearly enhanced with 
respect to the Bq = case. For different Bq, different peaks are enhanced. D{E) has values 
typically in the range ~ 10~^eV^^ to ~ 10^ eV^^ with enhancements up to ~ lO^eV"^ for 
energies E = — l.OeV to E = 1.0 eV. T{E) ranges between ~ 10^^'^ and ~ 10^^ for the 
same energy interval. The symmetry of D{E) and T{E) for positive and negative energies 
is apparent. The hopping parameters at the metal-torus contacts and in the device region 
were set to t-^Q^ = — 0.25eV and v = — S.leV. The Fermi level in metallic leads is chosen 
at -Efermi = 6.0 cV in rough accordance with physical values. 

Part (c) clearly demonstrates that improving the electronic coupling to the metallic con- 
tacts causes an increase in number of available states at the contacts into and out of which 
electrons can be scattered from and/or into the device region, as evidenced by a positive 
shift of the density-of-states D{E) with an increase in tj^Qp. %op niodified here from 
the default value %Qp = — 0.25eV to tj^Qp = — 0.5eV and t^^Qp = — 1.5eV respectively. 
Again, values for D{E) are in the range ~ 10"^ . . . lO^eV"^. Similarly, an overall symmet- 
ric increase in T{E) is observed in Figure 2d when increasing %Qp allowing close to unit 
T{E) ^ 1, for certain energies. 




(11) 



V vexp —A (r) • (rj — r^) 



(12) 
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B. Coherence in electronic transmission 



Up to now, results for D{E) and T{E) have been presented with the metalhc leads 
attached in a back-to-back configuration, i.e. the opening angle a is chosen as 180°. Having 
in mind analogous phenomena seen in fiat annular structures, the question as to whether 
the torus could act as a macro-molecular Aharonov-Bohm oscillator arises. When changing 
the opening angle a between the leads and applying a bias voltage across leads, the effective 
path length of electrons moving in the left or right branch of the nanotorus from metallic 
lead to lead will differ, opening up possible interference effects between the left and right 
electron wave functions. Due to the quasi-ballistic nature of transport in carbon nanotube- 
like structures, coherence in transmission through the nanotorus should be expected. The 
three-dimensional geometry of the nanotorus and its finite (non-zero) width will show a 



modulation in the oscillations of T(E) for different a observed in rings [2l|, |22, |23|, |2J, |26|, 



271 . l28l . I29I . ISOl . l3ll . I32I . l33l . l34j . l35l . l36l . 1371 . l38l . l39| due to differing possible electronic pathways 



with different major radii varying between R — a and R + a in the toroidal geometry. 

The recursive algorithm employed here easily allows for the modelling of a rotation of 
the position of one of the leads with respect to the other by determining the new positions 
of the four atomic contact positions for the rotated lead while the positions of the contact 
sites at the other fixed lead remain unchanged. A quick back-of-the-envelope estimate helps 
to determine the electron momenta k and thus for which de Broglie wavelengths A = ^ 
coherence in electronic transmission is visible before the finite size of the lattice constant of 
the underlying graphene lattice of a = lAA washes out any interference effects for higher 
momenta. The electronic energy E is roughly given hj E = Assuming E ^ O.OleV, 

calculate the occurence of interference minima in transmission in steps of Aa where de- 
structive interference occurs for pathlength differences As = (2n + 1) | with n = 0, 1, 2, . . . 
between the left and right toroidal branch: 

2m, 

R = 5.8nm M = 27ri? = 36.4nm I 



E ^ 0.01eV= ^ ^ A ^ 12.3nm n A/2 

^ • Aa^ 360° ■ = 60.6°. (13) 



The minimum angular difference at which the finite lattice size starts playing a role for 
a (3,3) armchair carbon nanotorus with 1800 atomic sites and for the dimensions used here 
is Aa = = 2.4°. Thus, for electronic energies E of ^ O.leV and larger, any coherence 
effects visible at lower energies will become negligible. 

The above predicted plateaus in T{E) due to constructive interference in transmission 
through the torus can be clearly deduced from Figures 3a and 3b. T{E) was calculated 
for E = O.OleV and E = 0.02 eV where plateaus according to Eq.(13) were predicted. 
The figures show several interference curves for an angular range of Aa = 45°, . . . , 315°. 
At larger energies, calculations show a gradual fiattening of the plateaus and eventually a 
disappearance of coherence effects at £^ = O.leV and E = 0.2 eV respectively. At smaller 
energies on the order of O.OOleV or less, T{E) appears essentially flat with small amplitudes 
~ 10~^ and less, quickly decreasing when turning on and increasing Bq. i?-field values from 
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T to larger fields of 5.0 T have been chosen to determine any dependence of the position of 
the plateaus as a function of flux threading the nanotorus plane. A substantial enhancement 
of T{E) by about one order of magnitude occurs for specific Bq values i?o ~ 1-0 . . . 1.5 T; 
the exact -B-field for maximum T{E) plateau shifts to larger values with increasing energy 
E. 

For comparison, the T{E) plots for non-resonant magnetic field values are shown enhanced 
in the insets. The resonant 5-field dependence of the transmission plateaus at the higher 
scale corresponds to the internal atomic scale of the torus, a = lAA. In this case, the 
magnetic phase ( = e**/*'' with $o = " creates a -B-field dependent electronic interference 
on top of the previously discussed interference effect due to different opening angles a. In 
simple terms, constructive and destructive interferences are generated by a magnetic phase 
change ( that accrue from lattice site to lattice site along the toroidal circumference to yield 
overall macroscopic phase changes for coherence in transmission from the left and right 
toroidal branches while at the larger scale this macroscopic phase change already occurs for 
electron hopping from lattice site to adjacent lattice site. Further details can be found in 
Subsection 3.3. 



C. IsD with magnetic flux; the torus as an Aharonov-Bohm oscillator 

Under a small voltage bias Vsd, Isd can be determined from Eq.(lO) by an integration 
over the transmission function T{E) convolved with the difference in the Fermi distributions 
of the two metallic leads. (Here the effect of the field on the leads has been neglected but 
should be considered in more comprehensive work.) Isd is plotted as a function of Vsd in 
Figures 4a and 4c for a number of different Bq while Figures 4b and 4d show the dependency 
of /5Z) on $ for two Vsd values. 

Linear volt age- current profiles for the small-bias case can be clearly shown in Figures 4a 
and 4c for back-to-back leads and 90° lead opening angle. The slope of the curves yields the 
conductance G xVsd with G as conductivity as predicted from Landauer-Biittiker theory in 



small-bias approximation [51l, l52|. The current can change by a factor of roughly 2.5 when 
choosing Bq = 0, 0.5, 1.0, 1.5, 2.0T. From this depiction no pattern in the dependency of Isd 
on Bq can be detected. 

Following the discussion of coherence patterns in transmission for different opening angles 
between leads, one could conceive of achieving similar interference effects by increasing 
magnetic flux $ = f Adr through the torus. The flux ratio ^ with respect to $0 can be 

c 

equated with a magnetic field ratio with Bq = 0.026 T. The current Isd is plotted here 
as a function of the flux ratio For both the 180° and 90° lead configurations, similar 
oscillations of Isd with changing magnetic flux are demonstrated. In each case, the strongly 
oscillatory curve appears to be composed of a superposition of two simpler oscillations in 
$ not shown here with $ as even or odd multiples of jq^o- Constructive and destructive 
interference in electronic transmission can thus be generated due to phase angle difference 
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with magnetic fluxes at the scale of j^. The flux period of these oscillations can be read 
off from the plots as A$ = ^$o- For comparison, magnetic flux dependencies are mainly 
computed and shown in the literature for persistent currents for two- and three-dimensional 
ring geometries [32, l33|, IS^, |43|, but the magnitude of Isd for example is comparable to the 
current magnitude expected in a two-dimensional graphene ring under similar conditions 
when re-scaling the results in 43|] to the torus dimensions here. Also, the general oscillatory 
behavior with changing magnetic flux appea rs to a gree with findirigs reported there and 



other places [22, l23|, M, M, M, M, M, M, M, \M, M, M, \M M, M 



To compute Isd (Vsd), T{E) is calculated symmetrically around the device Fermi level 
Ej between the energies E = — 0.2eV and E = 0.2eV for a step size of 5 x 10~^eV where Ej 
has been set to zero. T{E) is then integrated over this energy region with the difference in 
Fermi distributions according to Eq.(lO) to approximately yield Isd- For the small-bias case 
Vsd = 0, . . . , O.leV, contributions to the energy integration of Isd beyond the chosen energy 
range may be safely neglected due to the fast exponential drop-off of the Fermi functions. 



4. CONCLUSIONS 

The authors employ a recursive Green's function algorithm to calculate the electronic 
density-of- states D{E) and transmission functions T{E) of a thin toroidal carbon nanotube 
under a small bias voltage between laterally attached metallic leads. 

A strong enhancement of D{E) and T{E) is observed when the nanotorus is threaded with 
a magnetic flux $, particularly near van Hove peaks [46]. The enhancement is pronounced 
with increasing magnetic fleld. An increase of the hopping parameter t-^Q-p, which governs 
the effective electronic coupling at the metal-nanotorus atomic contact sites, also shifts D[E) 
and T{E) for all energies. 

Coherent electronic interference phenomena could be observed in T{E) as a function 
of the angle a between the metallic leads, and in the integrated source-drain current Isd 
as a function of magnetic flux. T{E) shows clear maximum and minimum plateaus for 
speciflc angular ranges where constructive or destructive interference of the electronic wave 
functions for electrons moving left or right in the torus from one lead to the other can be 
expected depending on the opening angle a. Due to flnite size effects in the nanotorus 
(i.e. the torus has a non-zero width) the constructive interference regions are broadened to 
plateaus over a range of angles and not simply peak-like maxima. Similarly, constructive 
and destructive electron interference appears in Isd when applying Bo where effective path 
length differences are induced by a magnetic phase angle e**/*°, different for the left and 
right toroidal branches. Isd measured in units of | shows strong oscillations in $/$o- The 
magnitude of Isd is comparable to the current expected in a two-dimensional graphene ring 



under similar conditions after re-scaling the results in [43[| to the toroidal dimensions. 

Natural extensions of the model and methods of this work are towards larger systems with 
more than tens of thousands of carbon atoms, and to include nearest and next-to-nearest 
neighbor contributions. For flat two-dimensional graphene nanoribbons, an opening of an 
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energy bandgap has been predicted [18|. Similar effects can be expected to appear in an 
expanded calculation for the three-dimensional geometries. Furthermore, a more realistic 
device setup with the metallic leads attached at the 'bottom' of the torus, i.e. the torus 
lieing fiat on the two planar leads will allow an increase in the number of metal-carbon 
contact sites and possibly lead to better and closer lead attachment. This could prove to 
be very interesting in the context of an applied microwave field with an additionally tilted 
static magnetic field due to the interplay between electronic excitations and persistent and 



bias-driven currents through the torus 40|, |4l|, |42] . It is conceivable that quantum control 



applications could follow when the induced dipole and solenoidal magnetic moments couple 
to electronic spin. 
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Figure captions 



Figure 1. Sketch of device setup with toroidal carbon nanotube and metaUic leads. 

Figure 2. Density-of- states D{E) and transmission function T{E) as a function of energy 
E: a. D{E) for different magnetic fields Bq. b. T{E) for different magnetic fields Bq. c. 
D{E) for different electronic coupling terms ij^op between torus and metallic leads, d. T{E) 
for different hopping terms t-^^^. 

Figure 3. Coherence in electronic transport. Plateaus in electronic transmission: a. Trans- 
mission function T{E) ai E = 0.01 eV as a function of magnetic field B for different angles 
a between metallic leads, b. Transmission function T[E) at E — 0.02 eV. 

Figure 4. Back-to-back leads: a. Source-drain current Isd as a function of source-drain 
voltage V5£i[eV] (small bias) for different magnetic fields Bq. Isd in units of |. Chemical 
potential at left/right lead: /ii 2 = Thermal energy: ksT = SOmeV. b. Source-drain 

current Isd as a function of applied magnetic field -Bo[T] {Vsd = 0.1 eV). 90° angle btw. 
leads: c. Source-drain current J^i^ versus source-drain voltage V5D[eV]. d. Source-drain 
current Isd versus magnetic field So[T]. 
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